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We study the implementation of global SO(4) ~ SU (2) l®SU (2)r symmetry in general potentials 
with N-Higgs-doublets in order to obtain models with custodial SO(3)c symmetry. We conclude 
that any implementation of the custodial SO (4) symmetry is equivalent, by a basis transformation, 
to a canonical one if SU(2)l is the gauge factor, U(1)y is embedded in SU(2)r and we require N 
copies of the doublet representation of SU(2)r. The invariance by SO(4) automatically leads to a 
CP-invariant potential and the basis of the canonical implementation of S'0(4) is aligned to a basis 
where CP symmetry acts in the standard fashion. We show different but equivalent implementations 
for the 2-Higgs-doublets model, including an implementation not previously considered. 
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I. INTRODUCTION 



Despite the enormous success of the Standard Model of Particle Physics (SM) in accounting for all the known 
experimental data on high energy physics to date, several theoretical and observational ingredients pressure the SM 
for extensions. One can list naturalness problems such as the hierarchy problem or observational facts such as the 
need for new sources of CP violation, the need for nonordinary weakly interacting stable matter (dark matter) or a 
t^-* ■ natural explanation for the smallness of nonzero neutrino masses. 

The simplicity of the scalar sector of the SM, however, has some unique features owing to its minimality. Because 
of the presence of one Higgs doublet, the Higgs potential possesses a global accidental 5*0(4) symmetry larger than 
the usual electroweak (EW) gauge symmetry SU(2)®U(1)y Even after electroweak symmetry breaking (EWSB), 
O r through the only symmetry breaking pattern, the SO(3)c portion of the global symmetry, usually called custodial 
JL ' symmetry 0-Q , remains in the potential. The symmetry is only approximate because it is respected by the gauge 
interactions only in the gy — > limit and by the Yukawa interactions only if the quarks in the same doublet have 
degenerate masses (it is then related to isospin symmetry). Nevertheless, such approximate symmetry guarantees the 
— '! gauge bosons W£ transform as a triplet of S0(3)c in the gy — » limit and explains the absence of large radiative 
. corrections to the p parameter from its tree- level value p^ — M^/M^ cos 2 9 W = 1, hence the name custodial 
' symmetry. 

, When we try to extend the SM, the most simple extensions of the SM consider some enlargement of the scalar 
I/"") ' sector of the theory since it is the least verified sector either because the only elementary scalar of SM, the physical 
Higgs boson, has escaped discovery so far or we have had no final confirmation about the mechanism behind EWSB 5] 
and the origin of all the masses of the elementary particles of the SM. The simplest way to extend the scalar sector 
j is to consider the replication of the Higgs doublet, leading to the N-Higgs-doublet models (NHDMs). The minimal 
version called 2-Higgs-doublets model (2HDM), was originally considered to implement the spontaneous CP violation 
mechanism Q , and more recently as simple extensions to address problems such as the little hierarchy problem Q , 
the presence and stability of dark matter candidates @, Q , and to provide new sources of CP violation when extended 
to a 3HDM 9[. This model has also been extensively studied in the last years as the effective scalar sector of the 
MSSM f]~ol - |l2i | . For example, Ref.[T3 presents a detailed study of the effective 2HDM potential that arises in the 
MSSM with explicit CP violation, after including radiative corrections. 

Another reason to study this type of model is the presence of a horizontal space formed by the TV-Higgs-doublets 
that possess the same gauge quantum numbers. An SU(N)h transformation in such space, then, only amounts to a 
reparametrization of the potential [i~3l - [i~7| . At the same time that the presence of a horizontal space makes the model 
more complex and less predictable due to the profusion of allowed couplings it also allows the study of the potential 
and its possible symmetry breaking patterns through the study of the orbit space of the gauge invariants [18l42l| . 

In the NHDM extensions, with N > 1, the global SO(A) symmetry does not arise naturally as an accidental 
symmetry but has to be imposed. Once it is imposed and leads to the custodial symmetry S0(3)c, it can play its 
original role to protect large radiative corrections to the p parameter |22H25l |. In the context of the 2HDM, the global 
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50(4) symmetry leads to the degeneracy of the charged scalar and one of the neutral scalars [2g|. In one version [23j . 
it is possible to construct a 2HDM model with custodial symmetry such that the state degenerate to the charged 
Higgs is a CP-even neutral Higgs, Mjjo = M H ± , instead of the more usual M^o = M H ± scenario. 

We intend to consider here the NHDM extensions of the SM possessing the custodial 5*0(4) symmetry which, in 
turn, can be broken down to the custodial 50(3)c symmetry that protects the p parameter. Larger global symmetries 
can be considered [27| but in general it will produce more pseudo-Goldstone bosons [28| in addition to the usual three 
would-be Goldstone bosons absorbed by the Higgs mechanism. We can separate our goals in this study into two 
categories: construction and identification. The first refers to the construction of general custodial invariant NHDM 
potentials and analysis of the models. The second goal refers to the methods of identifying the 50(4) symmetry if 
it is not manifest and concealed such as the CP symmetry acting in a general basis 0, E3, H3 • ^ is then important 
to study the interplay between the imposition of 5*0(4) symmetry and the reparametrization group SU(N)h- In 
particular, it is important to emphasize that the global 50(4) symmetry, like CP symmetry, is not contained in 
SU (2) l®U \\)y®SU (N) h as most of the symmetries usually considered for NHDMs are. See, for instance, the 
symmetries for 2HDM and 3HDM in Refs. 0, [2l|, [29|. The analysis of basis transformations is also crucial when 
extending some symmetries of the potential to the Yukawa interactions, since some symmetries of the potential might 
be broken by the Yukawa terms as well as the converse [22l l24l . [30L l3l|. 

In the context of the 2HDM, many studies [22, l24l. l25j| show that CP invariance and custodial symmetry are 
intimately connected. In fact, we will see that such connection can be extended to general custodial-invariant NHDMs. 
In certain cases the knowledge of the presence of the CP symmetry and how it acts on the fields, allows us to infer 
if 50(4) is also present as a symmetry. We can also borrow the various methods of identifying CP invariance 
independently of the basis [U, El, Compared to the SM, where global 50(4) symmetry, custodial 50(3)c 

symmetry and CP symmetry are automatic in the scalar sector, we will see that a NHDM potential invariant by 
50(4) symmetry and a potential invariant by CP symmetry differ only minimally in an easily identifiable manner 
for N — 2 and N — 3 Higgs doublets. For N > 4, a custodial-invariant potential is also more constrained than a 
CP-invariant potential but the difference is more involved. 

The paper is organized as follows: In Sec.|TT]we study one particular implementation of the global custodial 50(4) 
symmetry called canonical implementation and how it is related to the CP symmetry and basis transformations. In 
particular, in Sec. Ill Dl we show any implementation of the global symmetry 50(4) is equivalent to the canonical 
implementation. Section IIIII shows the various implementations for the 2HDM considered in the literature, shows 
their equivalence to the canonical implementation and also shows one different but equivalent implementation. We 
discuss some consequences of this work in Sec.HV] followed by the conclusions. To clarify the nomenclature regarding 
the various groups involved in the study of the custodial symmetry, e.g., 50(3), 50(4) and SU{2)l®SU{2)r 1 we 
refer to Appendix 1X1 where we list the group/subgroup structure of the global symmetries involved in the SM. Instead 
of distinguishing the various (custodial) symmetries by names, we will denote them by their group structure and use 
the isomorphic groups interchangeably, often considering their local structure only. We will often include a subscript 
a C" (custodial) to denote the groups, as in SO(A)c- 



II. CANONICAL IMPLEMENTATION OF SO(4) c IN NHDM 



Let us consider N Higgs-doublets of SU(2)l of hypercharge Y = 1: 4> a , a = 1,...,N. 
Let us also define a 2 x 2 complex matrix of fields [j, [j| 



a Ta 



<& a ^ = n , (i) 



where <j> a = e0* (e = £02) defines the tilde (~) operation on a doublet which transforms a Y = 1 doublet to a Y = — 1 
doublet. We can also define the (~) operation on any 2x2 complex matrix h as 

h = eh*J . (2) 
It is straightforward to see that such operation is trivial on $ a of Eq. (TT]), i.e., 

#a = $a • (3) 

We can define the Canonical Implementation (CI) of the custodial group 5J7(2)x,(g)5L T (2)_R by the universal action 
upon $ a , a = 1, . . ., N, 

$a -+ U L $ a Ul , U L G SU(2) L and U R G SU(2) R . (4) 
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The action Q sets the representation of $ a as (2,2) under the custodial group SU(2)l®SU(2)r. Therefore, each 
Higgs-doublet cj) a has the same transformation property under SU(2)l®SU(2)r as the SM Higgs-doublet. 

Within this section we will focus on the canonical implementation of the custodial group SO(4)c and any mention 
to the custodial group will denote this implementation. Within the 2-Higgs-doublet model, this implementation was 
named type I in an early work [22j . 

It is clear that the identity © preserves the action ^ since 

eU*e f = U, for any U G SU{2) . 
Using the operation ^fy we can write the previous identity as 

U = U, for any U G SU(2). (5) 

If the transformation (0} is promoted to a symmetry of the TV-Higgs-doublet potential, we say the potential is 
custodial-invariant. As in the SM, such symmetry is explicitly broken by the gauge interactions that are invariant 
only by the gauged SU(2)l®U(1)y subgroup. In the usual EWSB scenario, the three would-be Goldstone bosons 
that appear in the breaking of SU{2)l®SU{2)r — > SU(2)l+r coincide with the usual would-be Goldstone bosons 
of the breaking of SV "(2) l®U (l)y — > U(l) em . They are then absorbed as the longitudinal components of the gauge 
bosons W^^Z by the Higgs mechanism, such that no physical Goldstone bosons appear. Other symmetry-breaking 
patterns exist where physical pseudo- Goldstone bosons appear. That is the case of certain 2-Higgs-doublet potentials 
with global symmetries equal or larger than SO(A)c (see Sec. lII Bl and Ref.l27t). 

The correspondence between the representation (2, 2) of SU(2) l®SU(2)r and the vector representation 4 of SO(4)c 
guarantees we can construct a 4- vector out of the real and imaginary parts of the two components of <j> a . 

The usual correspondence is found by defining four real fields ?7 (a) = (77'°' ) = (?7Q a) , »7 (a) ) , rj {a> — (t]i ,T] 2 , ) i from 
the relation 

$ Q = „<•>! + = ( V % + S 1 + S) ■ (6) 

We can rewrite ([5]) in a more compact form 

*- = < ) e M . ( 7 ) 

by defining (e M ) = (1, iai, ia 2 , ia 3 ) . 

The custodial symmetry SU(2)l(3SU(2)r acting as @ induces in if a) the transformations 

<' -> R»v{U L , U R )n^ , R(U L , U R ) G SO(A) , (8) 

where the convention of summation of repeated indices is implicit for [a, v = 0, 1, 2, 3. It can be shown that R(Ul, Ur) 
defines a two-to-one mapping between SU(2)l®SU(2)r and SO(4)c with the explicit function 

R^(U L , U R ) = ±Tr[elU L e v U R ] . (9) 

Hence r/ <a) can be regarded as vectors in a four-dimensional euclidean space. We can also recover the SU (2) l+r 
SO (3) mapping for U R = U L ^. 

We can also define a 4-component complex vector (in a four-dimensional real vector space) that transforms more 
directly under SU(2) l ®SU(2)r. Let us define 



*"13J- (10) 

The action Q is then translated into 

Xa^U R ®U LX a- (11) 

The transformation property (fTTj) can be deduced if we rewrite (|4]) as § a ij — > ULii'U R jj,^ a i>ji [note that 

Xa = (e«)l2)($aii,$a2i,$ai2,$a22) T implies the transformation property \a -> {e<E)t 2 )(U R <E)U L )(^ <E)l 2 )Xa] and 
finally use the identity J3]). 

One can also write explicitly the linear relation between \a and 77 <a) as 

Xa = Brj^ , (12) 

where the transformation matrix B can be read off from ^ . With this change of basis it is possible to write explicitly 
how the generators of SU (2) l<E>SU (2) r acting on \a are related to the generators of SO(A)c acting on r/ (a) . Such 
generators can be found in appendix [Bl 
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A. Custodial invariants and explicit CP invariance 



For the construction of a renormalizable N-Higgs-doublet potential that is invariant by the custodial group SO(4)c 
in the canonical implementation we need to identify the minimal custodial invariants. 
The terms that are invariant by (j4]), up to mass order 4, are 

Tr[$t$ b ] = fy b + h.c, 

+ h.c. (13) 

It was necessary to use some reordering identities such as <j> d <fib4>t4>d = 4 , \,4 > d < pl4'c ~ ( t ) a4'c4'\ < t'd by making use of the 
identity eyejy = 5ikb~ji — o~uo"jk ■ There are no other invariants of the same orders and apparently independent invariants 
are reducible to the ones in Eq. (fT3"]) . For instance, det $ a = ^Tr[$t $ a ], and euiejji^ a ij^bi'j' = |Tr [$£$;,], due to 
©■ 

We immediately notice that the custodial invariants (fT3|) are invariant by the canonical CP (CCP) transformations 

CCP 

(j> a (x ,x) -> <l>Z(xo,-x), a=l,...,N. (14) 



Such transformations act linearly upon Xa (1101) as 



and upon rj i ' L) ([6]) as 



ccp i 6 

Xa "4 I _ £ q ) Xa, (15) 



. ^ . . ^ ) ^ (v ( o } , -V[ a> , V ( 2 } , -V^ ) ■ (16) 



We realize that 7] ( q } ,7]^ are CP-even fields while rf^^rj 1 ^ are CP-odd fields. Such transformation is in accordance 
with $ a -> $* = e$ a et. 

Since any custodial- invariant potential can be written as functions of the basic custodial invariants (|13l) . we can 
conclude that a N-Higgs-doublet potential which is invariant by the custodial group SO{A)c in its canonical form 
is automatically explicitly CP-invariant in its canonical form (| 14|) . This result was established in Ref.l22l for the 
2HDM potential considering the canonical implementation (type I) as well as a different implementation named type 
II. The latter, in fact, can be shown to be equivalent to the canonical implementation in a different basis 24, 25]. We 
will discuss basis transformations in Sec. lII CI and show in Sec. lIIDl that any implementation of SO(i)c in a NHDM 
potential is equivalent to the canonical implementation in some basis. 

We can ask the converse: Is an explicitly CP invariant NHDM potential automatically invariant by the custodial 
group 50(4)c? We know from the 2HDM potential [22|, [24], [25| that the answer is negative. The immediate following 
question is: what additional feature makes a CP-invariant potential also custodial-invariant? We will show in the fol- 
lowing, within the canonical implementation, that the additional condition is minimal for N < 3. The characterization 
for N > 4 will be seen to be more involved. 

For the 2HDM potential, the relation between custodial invariance and CP symmetry was already explored in 
Refs. [24], [Hj]. When basis transformations are allowed, such a relation can be related to the existence of a basis 
where the implementation of both the custodial group SO(4)c and the CP symmetry are the canonical ones. We will 
generalize this notion to general N-Higgs-doublets. 

To distinguish between a NHDM potential invariant by the CCP (IT4)) and by the custodial group SO(A)c (0]) we 
need to analyze the terms allowed by each symmetry to appear in a general potential. 

We know the most general renormalizable NHDM potential can be written in terms of the linear and quadratic 
combinations of the minimal gauge invariants [13] 



'6 • 



a,b = l,...,N. (17) 



Only AN — 4 of them are algebraically independent [l9| when the fields are considered as c-numbers. This fact affects 
the minimization procedure. 

Instead of working with complex quantities, we can consider the equivalent set of real quantities 

Re(^> 6 ), a<b = l,...,N, 

MKfo), a<b = 2,...,N. (18) 
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If we succeed in writing the gauge invariants (|18l) in terms of the custodial invariants () 13[) . then we can write a 
general custodial invariant potential. 

For the terms of mass order 2, we can check the correspondence is 

Re($&) = iTr[$t$ fo]j a <b = 1, . . . ,N, (19) 
= §Tr[$ t$ fo(7 3] , a <6 = 2,...,iV. (20) 

The term (I20p . however, is not custodial-invariant but transforms as the third component of a vector z (, = (z^) of 
SJ7(2) fl defined as 

4 = |Tr[$t$ 6CTA ] , A = 1,2,3. (21) 

In fact, the components 1 and 2 of z ab = are not U(l)y invariants as 

z\ h = ImttZefa) , z 2 ab = Re(^e^) . (22) 

Notice z aa = and there is no such term in the SM (N — 1). 

Both custodial invariance and CCP symmetry forbid any terms of the form below in the potential: 

Im(0t0 b ), a<b = 2,...,N, (23) 
Re(^ 6 )Im(<^rf) , a < b = 1, . . . , N, c < d = 2, . . . , N, (24) 

On the other hand, both custodial invariance and CCP symmetry allow the terms 

M<f>lM, a<b = l,...,N, (25) 
Re(^ 6 )Re(^ d ) , a < b = 1, . . . , N, c < d = 1, . . . , JV. (26) 



The only missing fourth-order terms are CCP-invariant but not custodial- invariant. They are 

lm{4,\cp b )lm{4>l(t>d) , a<b = 2,...,N, c<d = 2,...,N. (27) 
However, if combined with other terms as 

z ab -z cd = Im(<f>l(l> b )Im(<f>l<f> d ) + Re(-<#<M&<k + ^c^d) , (28) 
it becomes custodial-invariant. This term is contained in the second expression of (|13p since 

iTr[($t$ 6 )($t$ d) ] = Re( 0t b )Re(0t 0d ) _ Zab . Zcd . (29) 

Considering the properties z ba = —z ab and z aa = 0, we note that if some of the indices of (|28p are equal, excluding 
a = b or c = d, the term Im^^f^Im^J^) is effectively canceled out. For example, for b = c, 

Zah-zbd = Re(^0 b )Re(^^ d ) - (^ 6 Re($>d) ■ (30) 

Hence terms z ab -z c< i with two equal indices do not generate new genuine fourth-order custodial invariants besides the 
product of two quadratic invariants. 

The property (|3U)l means that custodial-invariant potentials with N = 2 and TV = 3 Higgs doublets do not contain 
terms of the form Im(0j l (/)f ) )lm(0j0 [ ;). The custodial-invariant 2HDM in its canonical form is, in fact, equivalent to a 
CP-invariant (real) potential without the quartic term Im 2 ((f>\<p2) [24l l25l |27T|. 

Therefore, we can state that a NHDM potential, with N < 3, invariant by the custodial group SO(4)c in its 
canonical implementation (j4]) is equal to a potential invariant by CCP (|14[) with the additional elimination of terms 
of the form Im(</>t cfr^Im^ltfid) ■ From a CP-invariant potential in its real basis, the latter elimination corresponds 
to the elimination of one term for N = 2 and 6 terms for N = 3. In this case, it is straightforward to construct a 
general 50(4)(7-invariant potential: only consider the terms in Eqs. (|23|) and (|2"6"|) . Such potential will have 9 and 27 
real parameters for the 2HDM and 3HDM, respectively. [A total of 2 + 2(r + q) + (r + q)(r + q + l)/2 real parameters, 
where r = N-l and q = N(N - l)/2.] 

One important remark concerning the description above is that for the full quantum theory considering radiative 
corrections, since gauge interactions do not respect the full custodial group SO{A)c : the effective potential will almost 
inevitably contain terms of the form Im^^t^Im^J^) generated by finite radiative corrections. But since CP 
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symmetry is a good symmetry of the potential and gauge interactions, the terms (j2"3")l and (|24p can only be generated 
from other sectors that violate CP, such as Yukawa interactions. 

For N > 4, the terms Im(</>£<^;,)Ini(</>J</>d) might be present in the potential contained in new genuine fourth order 
custodial invariants such as 

zi2-z 34 = Im((/4</>2)Im(<^4) + Re(-</>]>4<403 + <t>\fo<j>l<t>4) ■ (31) 

We note, however, that the custodial invariants z a h-z C( i, with a =/= b =/= c =/= d, still contain quadratic combinations 
of second-order custodial invariants Re ((/>£(/>(,), and we conclude they are reducible in some way. We are obviously 
interested in identifying the irreducible piece. 

The irreducible fourth-order invariants can be constructed most easily by using the real components ^ . We begin 
by rewriting the second-order invariants 

i) w 'f = |Tr[$t = Re (^ b ) . (32) 

The only SO(4)c invariant that can be constructed out of four 4- vectors rj {a) , r) {b) , r] {c) , rj {d) which is also independent 
of (|52"|) is the pseudoscalar 

2 abcd — ttivapVn Vu Va Vp ) 

where the e^ va p is the totally antisymmetric tensor in four dimensions obeying £0123 = 1 ■ We could have also written 

^abed = det V <abcd> f° r r) labad) ee (y a) |77 <b) |?7 (c) |77 <d) ) as a 4 x 4 matrix. The SO (4) invariant (J33]) is the genuine custodial 
invariant that only appears for 4 or more vectors (doublets). If any two of the indices a, b, c, d are equal the invariant 
is null. 

We can rewrite in terms of Higgs-doublets or z a b (|2ip as 

$L = im(0t&)M0£&) + M$0rf)M^0 c ) 

+ Im(<^ c )Im(<^ 6 ) , (34) 

= tabbed + Zac-Zdb + Zad'Zbc ■ (35) 

We note that the imaginary parts of the gauge invariants \m(<j)\(j)b) only appear through the combination 



zab-z c d = (v (a) -V (e) M b) -ri (d) ) - (r, M 'V w W b) 'V w ) 

+ I ( atld- (36) 

The invariant I a ^j cd changes sign by exchange of any pair of indices among (abed) and hence by any odd permutation. 
An even permutation of (abed), on the other hand, leaves I a ^j cd invariant. Thus each term ijQj is invariant by an 
A4 C SO(A)h discrete symmetry in each horizontal subspace spanned by {</> a , <^b, <A=, <^d}- 

Another useful formula is 

det(Imif( 1234 )) = (/W 4 ) 2 , (37) 

where if( 1234 ) is the submatrix of K, given by K a b — (frtcfia [3 [HI: constructed from the intersection of the rows 
1,2,3,4 with columns 1,2,3,4. 

Notice that the doublets live in the space C 2 and then at most two doublets can be linearly independent. Instead, 
when we are transported to R 4 through Eq. Q, four 4- vectors can be linearly independent, as they should be. For 
example, if 03 = a<f>i + P<f>2 for (f>i , <f>2 linearly independent, then still {rj {3) , 7y (1) , rj (2) } is linearly independent if Ima =/= 
or Im/3 ^ 0. 

For completeness, we can rewrite the SU(2)r vectors (j2"Tj) as vectors of SO(3)c: 

z afc = fjSV* ~ ?7oV o) + V ia> xr? (b) . (38) 

The most general custodial invariant NHDM potential for N > 4 can be written as the most general linear com- 
bination of terms of Eqs. (|23|). ([2"r?f and ([3~4"| . For the invariants I^ cd , there is only one invariant for each set (abed) 
with all distinct indices. For example, for N — 4, we know there is only one invariant. For N > 5, there will be (^) 
invariants. 

To perform the minimization of the potential with N > 5, however, we need to know additionally how many 
invariants 1^ are functionally independent among the (*T\ invariants. 
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B. Isotropy groups and symmetry breaking 



Let us analyze here what are the possible symmetry groups that are left invariant when the -ZV-Higgs-doublets 
acquire nonzero vacuum expectation values (VEVs), required from successful EWSB, in a potential invariant by the 
custodial group SO(4)c- The VEVs are realized as the absolute minimum of the Higgs potential and their isotropy 
groups define the spontaneous symmetry-breaking pattern. The specific patterns depend on the parameters of the 
potential but one can study the maximal realizable symmetry-breaking patterns from the possible isotropy groups of 
general VEVs with respect to both SO{A) c (global) and SU(2) L ®U(l) Y (gauge). 

Let us begin with the custodial group and consider N Higgs-doublets 4> a , together with their corresponding natural 
forms $ a and rj {a) , upon which SU (2) l®SU (2) r and SO(A)c act naturally as in Eqs. (U) and © respectively. We 
want to study how the possible isotropy groups He change as N grows. We will assume that all (<j> a ) ^ 0. Otherwise, 
it will coincide with the VEVs for smaller TV. In any case, successful EWSB also requires J2a( ( t > t. < t ) a) = w£/2, where 
r,. - 246GeV. 

For N = 1 Higgs doublet <\>\ (SM), SO (4) e is an automatic global symmetry of the Higgs potential and the 
custodial symmetry SO(3)c always remains in the potential after EWSB breaking because we can rotate away any 
corresponding VEV by a global symmetry SU(2)l alone and obtain 

($i>=«il, (39) 

where \f2v x = 246GeV within the SM. The isotropy group is H c = SU(2) L+R ~ SO(3) C - The VEV above automat- 
ically conserves CP. In terms of the real field rj (1> it is also obvious that any VEV (?7 (1) ) = 0, 0, 0) would leave an 
SO(3)c group invariant. 

For N = 2 Higgs doublets 4>ii4>2, the global symmetry SO(4)c is not automatic and has to be imposed in the 
potential. But once it is imposed, CP is an automatic symmetry of the potential. If we consider that EWSB proceeds 
through the VEVs ^ 0, ($ 2 ) ^ 0, then we can rotate to the form (j39j) by a global SU(2) L transformation. 
We can still apply an SU(2)l + r transformation [Ur = Ul in Eq. Q] which leaves (<f>i) invariant but allows us to 
write 



($ 2 ) = v 2 e~ lip ^ . (40) 

We have parametrized (r) <2) ) = (v2 cos <f2, 0, 0, — 1>2 siny^) for convenience. If ip2 ^ 0,7r, the isotropy group left 
invariant by and ($ 2 ) in Eqs. j39]) and (gDJ) is H c = U(1)t 3R +t 3R ~ SO(2) c ■ 

For N > 3 Higgs doublets, the isotropy group can be reduced to the trivial He = {1} or some discrete group. If 
($ 2 ) 7^ 0, we can rotate ($1), (<J> 2 ) to the forms ((521) and (gDJ) through global transformations. For ($3) ^ 0, we 
can use the U(1)t 3R +t 3R freedom to write 

($ 3 > = (r,< 3 »)l+( ?7 < 3, )za2 + ( ? ?3 3) )^3. (41) 

If (t7 2 3) ) ^ 0, there is no continuous symmetry left and the isotropy group is He = {1} or some discrete group. For 
the nonzero VEVs ($ a ), a> A, all the possible components are meaningful because there is no continuous symmetry. 

In terms of the gauge group SU (2) l(E)U (l)y , the possible isotropy groups H for N — 1, 2, 3, . . ., are known. For 
N = 1, H = U(l) em and the conventional form for (4>i) is 



«i>0. (42) 



It is the equivalent of Eq. (|39|) in terms of doublets. 

For N = 2, we can preserve the form of (|42[) . The isotropy group H is trivial if 02 is nonzero for the general form 



v 2 > 0, a 2 > 0, (43) 



where 6 2 is a nonzero complex number. For nonzero 02, this VEV corresponds to charge breaking vacuum (CB-V). 
This VEV contrasts to Eq. (|4"0"|) whose doublet version is 

(<fe)=f 2 ( e ° 2 ) » «2>0. (44) 

It corresponds to a vacuum that spontaneously breaks CP (CPB-V) in a real potential. In a general 2HDM potential, 
the VEVs (|43|) and (|44]) correspond to physically very different situations. In a custodial invariant potential, however, 
they are in the same SO(A)c stratum and have the same SO(2)c isotropy group with respect to SO(4)c- 
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In fact, in a custodial-invariant potential, a CPB vacuum and a CB vacuum arc always degenerate at tree level. 
In other words, if a VEV of the form (|44j) is a minimum of the potential, then there are infinite continuously 
related VEVs of the form (|43|) corresponding to the same minimum value of the potential. The reason is that 
V(U L ($i)U R , U L {^ 2 )U R ) = V(($i), ($ 2 )) and we can relate the matrices ($„) in Eqs. J3S]) and gDJ) to the following 
matrices through an SU(2)l + r transformation: 

Ul{$i)U r = ($!>, (45) 
Ul{$2)U r = w 2 (cos^ 2 l - sin^ 2 i<x-n) , (46) 

where n = (sin 6 cos a, sin sin a, cos 0) and 

U L = U R = e -^ 3 «/2 e -i CT2 e/2 _ ^ 47 ) 

The VEV l|46p corresponds to the form (|4"5| if we associate 6 2 = cos <p% +i sin <p 2 cos 9 and a 2 = sin </j 2 sin 6* for a = 3w/2 
in Eq. (|4"7|) . 

Another way of looking into the degeneracy between a CPB vacuum and a CB vacuum is to analyze the basis 
invariant that signals CB[l8j] : (r ) 2 — |(r)| 2 > 0. Only the equality signals a neutral vacuum. The variables r M , 
fj, = 0, 1, . . . , iV 2 — 1 are linear combinations of <j>\<j>b (fT7|) that will be properly defined in Eq. (jBT)) . In the simplest 
case of 2HDM, the r M variables are 2r = |</>i| 2 + |0 2 | 2 , 2r 3 = |0i| 2 — |0 2 | 2 , n = Re(</40 2 ) and r 2 = Im(</40 2 ), which 
allow us to write 

(r ) 2 - |(r)| 2 = K0i)| 2 K0 2 )| 2 - [Re(0l0 2 )] 2 - [lm(0^ 2 )] 2 . (48) 

Notice that the first two terms are custodial invariants while the third is not (f2T)]). The latter, Im(0{</> 2 ) = (zi 2 ) 7^ 0, 
signals SCPV in the real basis. Therefore, if a potential minimum obeys (r ) 2 - |(r)| 2 = but lm<0^ 2 > + 0, one can 
use SU(2)r to rotate (zi 2 ) in (|2"Tj) to decrease the 3-component and then increase (ro) 2 — |(r)| 2 in a 50(4)c-mvariant 
manner. 

If we recall the identity ([28]) for (cd) = (ab) = (12), 

|2 _ I J,. |2U_|2 _ r> Q 2/ -t 



one can see that 



zi2| 2 = |0i|1<fer-Re^ 2 ), (49) 



iwr = o (50) 



signals a neutral and CP-invariant vacuum in a custodial-invariant manner. In this case, the first two terms and the 
third term in Eq. (|48|) are independently null. 

The generalization of Eq. ([5U| to NHDM is given by 



a<b 

This condition ensures the symmetry-breaking pattern SO(A)c — > SO(3)c- 

If the symmetry-breaking pattern SO(A)c — > SO(3)c is realized in a NHDM potential, the residual symmetry 
SO(3)c guarantees that all the VEVs of the doublets are aligned as 

($ a )=v„l, o = l,2,.. .,JV, (52) 

restricted by 

2E u ^"i- ( 53 ) 

a 

with vl ~ 246GeV from successful EWSB. The residual symmetry also guarantees that spontaneous CP violation is 
not possible. 

It is important to emphasize that the degeneracy between CB and CPB vacua happens at tree level. Radiative 
corrections coming from the gauge sector and the Yukawa sector may break the degeneracy in the effective potential 
due to interactions that are not SO(4)c symmetric. Such tree-level degeneracy is a consequence of the symmetry- 
breaking pattern SO(A)c — > SO(2)c which generates 5 Goldstone bosons. If the low-lying vacuum, after radiative 
corrections, is neutral and then CPB, three of them will be absorbed as the longitudinal components of the massive 
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gauge fields. The remaining two scalars will combine to form a charged pseudo-Goldstone boson [28[ that acquires a 
small mass due to radiative corrections [13] ■ The same is true for any NHDM potential with only two nonzero VEVs 
that can be reduced to (|39|) and (|40|) . We are not interested in such a scenario due to the absence of the custodial 
symmetry that protects the p parameter and the lightness of such charged pseudo-Goldstone boson. 

The other scenario, i.e., the 2HDM scenario with a low- lying CB vacuum (with radiative corrections) or a NHDM, 
with N > 3, with a third doublet with nonzero VEV obeying Eg. (|41l) (tree-level), has no U(l) em symmetry and 
contains one more would-be Goldstonc boson that gives the photon a mass through the Higgs mechanism. This 
possibility is thus phenomenologically excluded. 

The actual parameter space that avoids the possibility of realizing ip ^ 0,7r in Eq. (|4TI)) or (ryf) 7^ in Eq. (T4"Tj) 
for the global minimum can only be mapped by constructing and analyzing the potential. For a 2HDM custodial- 
invariant potential with additional Z 2 symmetry, 2 — > —(j>2, both 50(4)c — > SO(3)c and SO(4)c — > SO(2)g 
symmetry-breaking patterns are realizable and the explicit parameter ranges for them can be found in Refs. pa. l27|. 



C. Basis transformation 



Since the N Higgs-doublets have the same gauge quantum numbers, any unitary change of basis will maintain the 
form of the kinetic terms and gauge interactions for the doublets but may modify the form of the potential as well as 
the Yukawa interactions without changing the physical content of the theory. It amounts to a reparametrization of 
the potential and Yukawa coefficients. 

The most general set of unitary transformations that commutes with the gauge group SU(2)l<S>U(1)y is given by 
the basis transformations 



u H e su(N) H 



(54) 



We are already factoring out the global ?7(1) factor because it corresponds to J7(l)y transformations. 
We can equally write the basis transformation law (|54|) for $ a (JlJ as [32] 



(t^r^iOc) 

= Re(U%)$ b +Im(U%)$ 



where 



i(-4>a 



$ a (-«)03 



(55) 



(56) 



Notice <£> a = <I> Q under @. If U H is a complex transformation, it is clear that linear combinations of {4> a } in the 
horizontal space of doublets will correspond to linear combinations of {$ a } and {$a}, i.e., only the set {& a } is not 
sufficient. From the definition (|56[) we see that a complex basis transformation implicitly induces a transformation 
in SU(2)ji and it shows that SU(N)h is not independent of SU(2)r. The transformation law for r] ia) can be easily 
extracted from Eq. (I5"5"j) . 

As extensively discussed [l3L Il7l . [29l | , basis transformations might conceal the presence of symmetries by changing 
how a symmetry acts. That is the case of CP symmetry for example [IH EH, 53] ■ That is also the case of the canonical 
implementation (0} of the custodial group SO(4)c- Hence, it is important to study how the representation (Q| of the 
custodial group SO(4)c is viewed in other bases connected by ([54]) . 

To that end, it is useful if we use a representation space where all the groups in question, i.e., SO(4)c and SU(N)h, 
act linearly through one matricial representation. It is better to treat 5*0 (4) c as in (fTTj) . 

If we define a complex column vector of size 47V (it lives in a 4iV-dimensional real vector space) by 



X 



(57) 



where 



the canonical representation (|TT|) of SU(2)t,®SU(2)r translates into 



(58) 



(59) 
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Notice 4> in Eq. (1571) denotes the application of the tilde operation on each doublet of 4> in Eq. (|58p . 
If we define new Higgs-doublets 4>' related by a change of basis (|54"|) . 

<t> = Ux4>'> ^ = (60) 
we can see the canonical representation (|59[) changes into 

rt 



$M%)<^^> (%)(.$)■ <6i) 



We can see the representations of SU (2) l is not modified and we will omit such factor within this analysis. The 
change in the representation of SU(2)r can be completely specified if we know the new representation of its generators 

{T' kR Y 

TU = (^ ^-r) T ™ ( Ul u*) > * = 1, 2, 3. (62) 
Obviously TkR = \o~k ® Ijv for the canonical implementation and we obtain explicitly 

rp, 1 [0 W 



where 



2 \iU 

T 3R = 5^3®liV = T3fl 



f7 = C/^i . (64) 

As expected, the representation of T$r ~ is not modified since the basis transformations ([54")) preserve the gauge 
structure. Only the representation of the generators Tir,T 2 r of the coset space SU(2)r/U(1)y is modified. 

We can also note that the representation for SU(2)r will be modified only when U\ in Eq. (|60f is complex. In fact, 
if we restrict the basis transformation group SU(N)h to SO(N)h, Ui is orthogonal, U — In in Eq. (jMJ) and the 
(representation of) generators (|55)) are the canonical ones. The converse is also true: the only subgroup of SU(N)h 
which preserves the representation of SU(2)r in Eq. (|6"5)) is the group of matrices Ui obeying UjU\ = In, which is 
equivalent to — Ui, i.e., they form the real subgroup SO(N)h- This fact is also in accordance with Eq. (l60l) which 
can be rewritten as Xa — (Ui)abX'b ^ or ^ SO{N)h- Hence, there are infinitely many canonical implementations of 
the custodial group SO(4)c related by SO(N)h- 

We can interpret the previous fact as the realization that only the real SO(N)h subgroup of the basis transformation 
group SU(N)h commutes with the custodial group SO{A)c- This also means that the whole global symmetry of the 
potential can still be enlarged by discrete or continuous subgroups of SO(N)h without affecting the custodial group. 



D. Equivalence among different implementations 



We will show in this section the following: Any implementation of the custodial group SU (2) l^SU (2) r/ Z2 on a 
potential containing N Higgs-doublets corresponds to a canonical implementation (0]) in some basis provided that 

1. the Higgs kinetic term is invariant by the custodial group; 

2. the association T 3R = \Y fix U(l) Y C SU(2) R ; 

3. SU(2)r acts through a reducible representation of N copies of the fundamental 2-representation. 

We already analyzed the converse in Sec. lII Cl i.e., a canonical implementation of the custodial symmetry satisfies 
conditions 1, 2 and 3 in any basis. 

To prove the statement, we begin by considering x i n Eq. (|57[) as the appropriate representation space of the N 
Higgs-doublets <p a . It has the same degrees of freedom of considering the real and imaginary parts of each component 
of each doublet as independent (4iV real fields). Then we note that the representation of SU(2)l is fixed by the 
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gauge interactions and each doublet <f> a transforms according to the same fundamental representation 2. If Ul is one 
element of SU (2) l in the fundamental representation, it acts reducibly on X as 

X SU ^ )L {t2N®U L ) X . (65) 

The action above defines the explicit representation D(Ul) — ^2N®Ul- 

For the SU(2)r factor, we want the representation D(Ur) of Ur G SU(2)r to commute with any Ul G SU(2)l, 
i.e., [D{Ul), D(Ur)] — for any Ul, Ur. From the irreducibility of 2 of SU(2)l, the only possibility of the action of 
SU(2) R is 

X SU ^ R (D« N \U R )®1 2 ) X . (66) 

We need to specify the representation D ( (Ur) in the representation D(Ur) = D ( (Ur)(B)12- The following consis- 
tency condition should be additionally imposed on such representation in a later stage: 

(e®l N ®e)x* = X- (67) 

For the SU{2)l factor the condition above is automatic. Note that Eq. (|67| means that complex conjugation acts as 
a linear transformation on \- 

We begin now to impose condition (1), which requires the form-invariance of the Higgs kinetic term 

d^d»<P=±d^d>*x, (68) 

under the global symmetry SU(2)l®SU(2)r. The invariance under SU(2)l is automatic. The invariance under 
SU(2)r implies the representation D ' should be unitary, i.e., 

(lf\u R )j 'd (2N \Ur) = 1 2N , for all U R G SU(2) R . (69) 

Recall that any SU (2) representation is equivalent to a unitary representation because it is a compact Lie group. 

The representation D ( ' N) of the Lie group SU(2)r also induces a representation of the respective Lie algebra su(2)r 
through the exponential mapping. The unitarity of D ) for SU (2) then induces a Hermitian representation for su(2). 
We can then work on the group and the algebra interchangeably if we are only interested in the connected part of 
the group modulo the center. We know su(2) is a three-dimensional Lie algebra and the Pauli matrices T k = \ok, 
k = 1,2, 3, can be used as a basis of the su(2) abstract Lie algebra obeying [n,Tj] = itij^k- We can also regard the 
explicit form of {rfc} as the basis for su(2) in the fundamental representation. Using the basis {r^}, a general element 
in the su(2) real algebra can be written as r^afc, where a^, k = 1, 2, 3, are real coefficients. The exponential mapping 
then reads 

(2iV) (2W) 

Ur = exp(iT k a k ) , D (Ur) = cxp(L> (iT k )a k ) . (70) 

We want to embed U(l)y into SU(2)r by requiring condition (2), i.e., by associating hypercharge to the unique 

diagonal generator of SU(2) R : T^r = . The action of hypercharge is fixed by the gauge quantum number 
assignments 

x VQ}y e iHY x= ( e ^a 0ljv0l2 ) x , (71 ) 
Thus, within the iV-Higgs-doublet sector, we obtain 

T 3R = T 3 ®l N ®~k2 ■ (72) 

and 

(2N) 

D (r 3 ) = t 3 ®1jv • (73) 

Now we can identify U(l)y ~ U(1)t 3R - Notice we could have associated T 3 r — —Y/2 and this would only change 
Eq. ([731 to D 2N (ts) = — T3®ljv, and the rest of the generators of SU(2)r accordingly. 

We also know all the representations of SU(2)r are totally reducible and our requirement (3) of the presence of N 
irreducible 2-representations implies that there is a 2Nx2N nonsingular matrix U 1 for which 

U'-'d'^XUrP' = U R ®1 N , for all Ur g SU(2) r . (74) 
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In the algebra Eq. ([71)1 translates into 

U'- 1 D (2N \r k )U' = T k ®l N , for k = 1,2,3. (75) 

Notice the usual block diagonal form (direct sum) of N copies of irreducible 2-representations for r k is l]v<8>Tfc = 
diag(rfc, Tfc, . . . , Tfe). But the latter can be transformed into the form (|75|) by a similarity transformation. 
Now, by using Eqs. (|T3"]) and Eq. (1751) . we obtain the condition 

[[/',r 3 <8)l J v]=0. (76) 

Since T3<S)1n is already diagonal, U' has to have a block diagonal form, each block corresponding to distinct eigenvalues 

Of T 3 ®t N 

Ui 



u '={ou 2 )- ^ 

The matrices U\,U 2 are nonsingular complex matrices of size NxN. Additionally, the unitarity of the representation 
D l ) (|69[) or the hermiticity of ([75)1 implies U\, U 2 have to be unitary matrices. 
The consistency condition ((67)) yields 



(e®lN)(D l2N Xu R ))*{e®l N ) f = D <2N \U R ) . (78) 

By using (|77[) we can rewrite 



-t N \ fu; \ > fuj \ / 1 N \ 



^)(^)(? ^) 



which reduces to 

^ D7°i7 9 ) (c/ ^ lw) = 1 uju 2 ) ■ ( 8 °) 

The identity ([5]) was assumed. 

Since the previous identity has to be true for the algebra as well, we can use Ur —> ir± from which we conclude that 

UjU 2 = UjUi . (81) 

The same condition arises from Ur — > ir 2 . The i-factor in it\ , ir 2 should be present to preserve the condition (|67p , 
which should be valid in the group representation but the latter is related to the representation of the algebra by the 
relation ((70]) . 

The explicit representation for r k (|75p is then 

ifM = i ( » "f ) (82) 

where 

C/ = U 2 U\ . (83) 

Notice the matrices in Eq. (|82p represent appropriately the algebra [7}, ry] = ieijkTk- 
We can check U is symmetric by using (fBTj) 

C/ T = t/j'Kj = U?UjU 2 Ul = U . (84) 

If we compare Eq. (f8"2"|) to Eq. (|M|) , we conclude immediately that the matrices in Eq. ([521 correspond to the generators 
of SU(2)r implemented in some basis cf>' = U\<f>. The unitary basis transformation matrix Ui should satisfy 

UjU 4 = U , (85) 

where U is given by (153")) . It was shown in Ref.[l4| (Appendix B) that the decomposition (|551 is always possible for 
a unitary symmetric matrix U. Thus, there is always a basis for which any implementation of the custodial group 
obeying the requirements 1, 2 and 3 can be cast into the canonical form ([I"]), 
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E. Custodial-invariant potential in the canonical form 



In this section we want to analyze the most general NHDM potential invariant by the custodial group SO(4)c in 
the canonical implementation (CI) and seek general features. Since any implementation of SO(4)c is equivalent to 
the CI in some basis, one can always use the latter basis to treat the potential. There are infinitely many such bases. 

We begin by writing the most general renormalizable NHDM as [lj, [l8| 

V = M M r M + A^rv , (86) 

where r M are variables that depend bilinearly on the fields (doublets) and {M^,A^} are coefficients. The convention 
of summation of repeated indices is implicit with Euclidean metric instead of the Minkowski metric used in Ref . fl7l. 
In this way, all the indices referring to the variables r M will be written as lower indices. 
The variables r u are defined as [H El 

r M = (Tjabtlfo , (x = 0,l,...,d. (87) 

The NxN matrices T M are 



where {A^} arc the d = N 2 — 1 hermitian generators of SU(N)h in the fundamental representation, obeying the 
normalization TrfAjAj] = 25ij. 

We can choose {Ti} to be the generalization of the Gell-Mann matrices of SU(3) to SU(N). We separate them into 
three classes: diagonal matrices {hi}, off-diagonal symmetric matrices {S a b}, and off-diagonal antisymmetric matrices 
{Aab}, with r, q, q elements, respectively [r = N — 1, q = ^N(N — 1)]. The matrices {hi} correspond to a basis for 
the Cartan subalgebra while the S a b and A a b correspond to the symmetric and antisymmetric combinations of ladder 
operators labeled by the pair (ab), a < b, a, b = 1, . . . , N. More explicitly, we can write 

Sab = \{e a b + e ba ) and A ab = j- t {e a b - £ba) , (89) 

where e ab are the canonical matrices obeying (e a b)ij — 5 a iSbj- We can follow, for instance, the ordering 

(ab) = (12), (23),..., (N -1,7V), (13), 

...,(N-2,N),...(1,N). (90) 

For example, for SU(3) we would have {h\,h 2 } = {A3/2,As/2}, {S12, S23, £13} = {Ai/2, Ae/2, A4/2} and 
{A12, A23, A13} = {A 2 /2, A7/2, A5/2}. The ordering (|M| follows the ordering of the lower to higher heights of the 
roots and the first r of them correspond to the simple roots. 

We will denote the variables corresponding to {hi, S a b, A ab } respectively by 

{*"(<) .^(aft),^}. (91) 

Therefore {^m} are combinations of 4> a <j> a while 

r {ab) = Re(cf>l<f> b ) , = lm(0^ b ) , a < b . (92) 

Then the elimination of the terms (|2"3")l and (|2~1T) corresponds in Eq. (|5d]1 to 

= for fi = {ab) (93) 



and 



= for fj, = (ab) or v — (ab) , (94) 



where or means one of them exclusively. The restriction (|94p leads to the block diagonal form for the dx d submatrix 
A = (Aij): 

A= (^ pxp ® pxq ) , (95) 
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where p = r + q = N(N + l)/2 and O pX g denotes the null matrix of the shown dimension. 
For N < 3, we have to impose additionally 

A Mi , =0 for /i = (ab) and v — (ab) . (96) 

The restriction (|96| sets B qxq — qxq in Eq. ((95|) . The most general custodial invariant 2HDM and 3HDM potentials 
are then given by flSBJ) with the restrictions p?]). (jMj) and CT : 

F= A Vm + A^r M r„, (97) 

;u=0,(i),(a6) /n,i/=0,(i),(a6) 

where the summation is explicitly shown. 

In this case, we can confirm the direct connection between the presence of only one charged would-be Goldstone 
boson and a CP-conserving vacuum. Such connection follows from the complete absence of the terms Im(^^(,) in the 
potential (j97|) . from which we immediately conclude that independently of r M , 

dV 

= 0, a<b = l,...,N. (98) 



Therefore, we can immediately write the extremum equations as [if 



dV x , dV 

— - = (M) o6 ^ i = ^ ^-(r^Bi^M. (99) 

0al ^0,(i),(cd) ^ 

We know (M), the matrix M computed for the minimal values (VEVs), is the squared-mass matrix for the charged 
scalars for a neutral vacuum [18|. including the charged Goldstone bosons. If there is only one (would-be) charged 
Goldstone boson, then the VEVs corresponding to the second components of each doublet, (w a ) = {4>a2): form the 
unique eigenvector of (M) associated with the null eigenvalue. We can see from Eq. (|99l) that M is real and symmetric 
because the sum only contains generators T M that are real and symmetric. Since the eigenvalues of (M) have to be 
real, the eigenvector (w) is also real and then spontaneous CP violation is excluded [33| . Obviously, this happens 
because we are implicitly selecting the symmetry-breaking pattern SO(4)c —> SO(3)c- The direct connection above, 
however, is not apparent for custodial-invariant NHDM potentials with N > 4. 

The independence of the potential on any r,^ (j9"5j) also explains why CB and SCPV vacua are degenerate for the 
custodial-invariant 2HDM. In this case, the potential does not depend on r^, which can be changed continuously and 
independently in that direction (within the surface and interior of the lightcone in Ref. [2lj ) without changing the 
potential value. The physical range of the variables r M ((57)) in the 2HDM corresponds to all the points inside and on 
the surface of a cone (lightcone) defined by r\ — r 2 = 0, where r 2 = r\ + r| + r\. The VEVs lying on the surface of 
such a cone correspond to neutral vacua while the interior points correspond to charge-breaking vacua. Then, if there 
is a neutral but CP-violating vacuum (r^) obeying (rg — r 2 ) = 0, with {r-z) ^ 0, the vacuum with the same values of 
r o, r i, r 3 but (ra) = lies at the same depth of the potential but it is charge breaking since (r 2 — r 2 ) > 0. The same 
conclusion can be reached for 3HDM if we can find independent directions to move in the subspace of r2 , r$ , rj in the 
orbit space [l9| . 

For N > 4, the most general custodial invariant NHDM potentials are obtained by adding to (|9"T|) the term 

AV^ ]T Xabcdlitd- (100) 

a<b<c<d 

From the definition of 7 (4) (34]) we can establish the correspondence 

A Q 6cd = 2A^^g) = -2A(_ = 2A^^ , (101) 

for a < b < c < d. Equation (|10ip restricts A^^ for fi,v = (ab). For N = 4, the expression (|100[) contains only one 
term Z^L • 

Notice that for N > 4, Eq. (IM|) is not valid. Because of Eqs. (jlOOp and (|101D . however, the coefficients A^ in the 
sector of T^Xr^ is a sparse matrix. For example, for N — 4, in the subspace of {^(12)? r (T3): r (Ti) ; r (23): r (24)' r (34)}j 
the submatrix B qxq in Eq. (|95p can be written as 



( A (ab){cd] 



#6x6 - (A/r^z-ft — A1234 



_ O3X3 B\ 

2)) ^ 1234 ^ ^ 03><3 
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B l = 0-10. (102) 



It is clear that the block matrix in Eq. f| 102[) has eigenvalues ±Ai234, each one with multiplicity 3. The explicit 
construction for N — 5, reveals that the corresponding block matrix has eigenvalues 0, A, — A, with multiplicities 
4, 3, 3, respectively, where A is a positive number that depends on the 5 parameters \ a bcd, a < b < c < d. 



III. DIFFERENT IMPLEMENTATIONS IN 2HDM 



We analyze in the following the different, but equivalent, implementations of SO (4) c in the 2HDM potential. The 
first two implementations, Secs. lIII Al and MI Bl were already considered in the literature but we show in Sec. lIII Cl 
that they can be generalized. 



A. Twisted implementation 

It can immediately be seen that the "twisted" implementation of the custodial symmetry in Ref.[23| is just the 
canonical symmetry Q imposed in a basis 

#, = e w »0 o> a =1,2. (103) 
It corresponds to a change of basis in Eq. (|60j) with 

ttf =diag(e"-), (104) 

where we can choose ^2 a a = from U(1)y invariance. 

This implementation can be straightforwardly generalized to iV-Higgs-doublets with general N > 2. 

For the particular case of 2HDM, two apparently inequivalent implementations of the custodial group were presented 
in Ref-HH- The twisted version was shown to allow for the degeneracy between the charged Higgs -ff ± and the CP-even 
scalar H°, Mh = Mg±. However, as shown in Sec. Ill Dl any implementation of the custodial group is equivalent 
to the canonical one in some basis. In the specific case of the twisted version of Ref.[23|, the required change of 
basis is the rephasing transformation of (I103p . After the appropriate basis change, the canonical CP transformations 
automatically constitute a symmetry of the Higgs potential. Let us call such automatic symmetry CP1. Thus the 
"twisted" implementation corresponds to an imposition of an additional CP2 symmetry, orthogonal to CP1 in the 
sense of Ref.[t7|, apart from the automatic CP1 symmetry. With respect to the CP1 symmetry, the CP2-even scalar 
H° is actually CPl-odd. Obviously, one can choose one of the CP symmetries to be a true symmetry of the theory 
by extending it to the Yukawa sector. Only then, is it possible to define (approximately) the scalar degenerate to 
as CP-even, the CP symmetry being the CP2 transformation. A similar comment has been already made in Ref.[24l. 



B. Type II implementation 

For the specific case of the 2HDM, various apparently different implementations of the custodial group can be 
devised. An implementation called type II [22| follows from the action of SU{2)l®SU{2)^ upon 

M l2 ={^x), (105) 

instead of the action on $i , $2. This implementation was first considered in Ref.[22| as distinct from the canonical one, 
which was named type I. However, Refs. [24j, |25j show that they are in fact equivalent to a canonical implementation 
in a different basis. 

To complement such discussion, let us analyze the equivalence between type II and canonical implementation from 
a slightly different perspective. 

The matrix M12 is a general complex 2x2 matrix, different from $ a which obeys ©. In terms of independent 
degrees of freedom (real fields), $1 and $2 together has 8 degrees of freedom, the same number as for A/12. The 
action of SU(2)l®SU(2)r on M 12 in Eq. (|105[) is, however, reducible because of the property ([5]). The irreducible 
pieces of M12 can be separated, respectively, as the even and odd parts under the tilde operation @: 
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M «-, _ (106) 

We can also think the tilde operation as a linear transformation in the space of Xa (|10j) with two eigenvalues, ±1, that 
separates the irreducible representations of M\i. 

In other words, we can define doublets 4>[ and 2 associated, respectively, to 

$ 2 = -iM& = . (107) 

or, equivalently, 

1 

72 



M 12 = —(a^ + ib^e^, (108) 



where a M and b^ are 4- vectors similar to ?y (a) . The explicit relation between the doublets <p a in (I105p and <p' a is given 
by 



V2 

h - 02 



V2i 



(109) 



One can see, as expected, that the type II implementation on (|105[) corresponds to the canonical implementation 
on cj)' a through &' a constructed as Eq. ([1]). It should be remarked that the basis change in Eq. (I109p connecting type 
II implementation and canonical implementation is not unique because any additional SO{2)h transformation on <p' a 
still preserves the canonical implementation of SO(4)c- 



C. Other implementations 

We have shown in Sec. Ill Dl that all implementations of the custodial group are equivalent to the canonical im- 
plementation in some basis. For the 2HDM potential we can show another practical way of implementing custodial 
symmetry, in clear analogy with type II implementation, without resorting to basis change. 

In this class of implementation we define 

My i = (0i 1 01 ) , (110) 

where 

0i = Cl 0i+ C2 2 , \ Cl \ 2 + \c 2 \ 2 = 1. (Ill) 

We can easily see (ci,c 2 ) = (1,0) [in this case $ 2 will be also necessary] and (ci,c 2 ) = (0,1) correspond to the 
canonical and type II implementations respectively. 
The action of SU(2) L ®SU(2) R is given by 

M V i^U L M vl U ] R . (112) 

By analyzing the infinitesimal transformation of SU{2)r on 0i and 2 , we can conclude that such implementation 
is equivalent to the canonical implementation on (</>' 1; <f>' 2 ) given by 




(113) 

One can check the previous transformation matrix is unitary and symmetric. 
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IV. CONSEQUENCES AND CONCLUSIONS 

We have shown in Sec. HID I that any implementation of the custodial SO(4)c ~ SU(2)l®SU(2)r group in a 
NHDM potential is equivalent to the canonical implementation in some basis if SU(2)l is the gauge factor, f/(l)y 
is embedded in SU(2)r, and we require N copies of the doublet representation of SU(2)r. On the other hand, we 
have seen in Sec. Ill Al that a potential invariant by SO(A)c in the canonical implementation (CI) is automatically 
CP-invariant in its canonical form (CCP). Therefore, any NHDM potential invariant by a global SO(A)c symmetry 
is also automatically CP-symmetric independently of the implementation. Moreover, the basis for which SO(A)c 
implementation is canonical is automatically aligned to the basis where the CP transformation on the doublets is 
canonical. If one requires the custodial symmetry SO(3)c to remain after EWSB, spontaneous CP violation is 
automatically excluded as well. 

One can then take advantage of such result to use explicit CP invariance as a necessary criterion to identify global 
SO(4)c invariance in a NHDM potential, even if its manifestation is not explicit. Several tools based on basis invariant 
conditions were already developed for identifying CP-invariance in a general NHDM potential [l3|, [3 [l?} • One can 
also make use of the fact that a CP-invariant NHDM potential can be always written in the so-called real basis (or 
the canonical CP basis [T3|) where all the parameters of the potential written in terms of the gauge invariants (f>\^)b 
are real[T3|. In other words, there is always a basis where the CP symmetry manifests canonicallv (fT4")l . 

To find such basis for a general CP-invariant NHDM potential with N > 3 is a difficult technical and unsolved 
problem [lTj]. For the 2HDM, the change of basis from a general basis to the real basis can be explicitly given [17j . For 
the 3HDM, necessary and sufficient criteria for CP invariance can be devised but the systematic method to find the 
real basis is lacking |17| . In both cases, once we write the CP-invariant potential in its real form, then identifying if 
the potential is <S , 0(4)c , -symmetric is a straightforward task because any real basis is also a basis where the action of 
the SO(A)c symmetry is canonical. Thus, it suffices to check if the potential is entirely written in terms of Re(<%<^,) 
in its linear and quadratic combinations, just as in Eq. (1971) . Such a criterion is invariant by the residual SO(3)h 
reparametrization freedom that remains. 

For the specific case of the 2HDM, basis-invariant conditions to test SO(4)c can be written 25]. If we follow the 
methods of Ref.ll7l we can consider the parameters {M = (Mi), Aq = (Aot),A = (Ay)} as two vectors and one 
rank-2 tensor of S0(3)h, respectively. These parameters were defined in Eq. (|56"1) . CP transformations then act 
through ordinary reflection in three dimensions along a direction k-cp called CP-rcflcction direction, which might not 
be unique. CP invariance then is equivalent to the existence of kcp such that 

kcp-M = , kcp-Ao = , kcp is an eigenvector of A . (114) 

In addition, the potential is SO (4) c -invariant if, and only if, 

kcp is an eigenvector of A associated to the eigenvalue . (115) 

Obviously det A should be null [25| . The possible directions for kcp were given in Ref.[l7l but when M, Ao are not 
null and nonparallel vectors, kcp = MxAo is certainly one CP-reflection direction. 

For 3HDMs, a systematic procedure to find the real basis was not completed but we can still write the necessary and 
sufficient conditions to have SO(A)c symmetry. If the potential is CP-symmetric, there should be a three-dimensional 
CP-odd subspace t q of the adjoint space, the 8-dimensional euclidean space where r = (fj) lives, such that [l?} 

M _L t' q , Ao -L l q , is an invariant subspace of A. (116) 

In the real basis, t' q is spanned by {A2, A5, A7}. The potential is additionally S , 0(4)c , -invariant if, and only if, 

tl is contained in the nullspace of A . (117) 

For TV > 4, a real basis may not coincide with a basis of CI of S0(4)c- More specifically, CCP invariance (real 
basis) requires the parameters {M^A^} to obey some but not all of the restrictions for SO(4)c invariance in the 
CI. The restrictions for CCP invariance coincide with Eqs. (|93p . and (|95[) . but S0(4)c invariance additionally 
requires that B qxq in Eq. (|9"5j) should be compatible with the structure in Eq. (|100D . 

For example, for N = 4, a CP-invariant potential can be brought to the real basis where we can identify the block 
matrix B qxq in (A Ml/ ). If the potential is also S , 0(4)c-invariant, the structure of B qxq should be of the form in 
Eq. (|102p (CI) or S0(4)h equivalent. The group 50 (4) # appears as the real subgroup of the reparametrization group 
SU(A)h which preserves the representations of 50(4)c and in fact it is the only compatible reparametrization group. 
Now, an element g of the group 50(4)# acts on B ex6 as 

#6x6 -> D 6 (g)B 6x6 D 6T (g), g G SO(A) H , 
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where D 6 is the 6-dimensional adjoint representation of g. Certainly, from Eq. (|102[) . 50(4)c invariance requires 
-86x6 to have only two eigenvalues ±A with multiplicity 3 each. The problem is that such condition is not sufficient to 
guarantee that there is a matrix D 6 (g) that transform Bq x q into the form (|102[) because a 5*0(6) orthogonal matrix 
is necessary in general. Therefore, additional conditions are necessary. 

In the same way, we can extract the necessary condition for N = 5: B qxq in the real basis should have eigenvalues 
0, A, — A with multiplicities 4,3,3 respectively. Analogous necessary conditions for general N > 4 could be found by 
studying the terms in (jlOOj) . 

As for the general problem of constructing the most general NHDM potential invariant by the global 50 (4) c 
symmetry, we managed to solve it in its simplest form: we can write the most general custodial-invariant potential 
before EWSB. The details of EWSB and the parameter range to achieve the desired symmetry-breaking pattern 
SO(4)c — > SO(3)c will be left for further work. For the ^-symmetric 2HDM, the parameter range is known since 
Refs.j2H3|. 

In any case, we confirmed that the phenomenologically interesting symmetry-breaking pattern is SO(A)c ~ > 50(3)c- 
In this case, the custodial symmetry 50(3)c ~ SU(3)l+r remains as an approximate symmetry at the electroweak 
scale and protects the p parameter of acquiring large radiative corrections quadratic in the masses of the scalars. 

The other viable 50(4)c — > SO(2)c symmetry-breaking pattern does not have the custodial symmetry and its 
viability in realistic models can not be treated at tree level because the CP-violating but neutral vacua are degenerate 
with the charge-breaking vacua. Only if radiative corrections are taken into account, can we check which models 
have one neutral vacuum lying deeper than the charge-breaking vacua. This fact was not sufficiently emphasized in 
previous works [24l [2I \2l\ . 

Because the symmetry-breaking 50(4)c — > SO(2)c is possible and spontaneously CP-violating and charge-breaking 
vacua might be degenerate, we can conclude that the imposition of the custodial group SO(4)c before EWSB and the 
requirement of SO(3)c custodial symmetry are separate conditions in NHDMs (with N > 1). The imposition of the 
global SO (4) c symmetry certainly implies explicit CP invariance in the NHDM potential but the possibility of having 
spontaneous CP violation is still present with the possibility of the pattern 50(4)c — > SO(2)c- 

An interesting question that arises through the study of the global 50 (4) c* in NHDMs is the role played by global 
symmetries in the potential. For the 2HDM, with Z2 symmetry, different global symmetries ranging from the usual 
Gsm to 50(8) were studied in Ref.l27l In general, for a global symmetry larger than 50(4)c, (pseudo) Goldstone 
bosons appear [28|. For example, the global symmetry group Gsm <8> SU(2)h is broken down to U(l) ® U(l) em 
generating three would-be Goldstone bosons and two true Goldstone bosons associated with the breaking of SU{2)h- 
These true Goldstone bosons are neutral. Reference [27| also explored the possibilities of having 50(4), 50 (4)® 5*0 (4) 
and SO (8) as global symmetries. The case of the 50(4) group is the usual custodial symmetry extended to 2HDM. 
The case of 50(4)(g)50(4) is the unlocked case where the custodial symmetry acts independently for both doublets. 
In these cases the global symmetry is not an independent factor but contains the gauge symmetries of the SM. A 
detailed study of possible approximate symmetries in NHDMs is performed in Ref.(3J. 

The maximal continuous global symmetry in 2HDMs is SO (8) contained in 0(8). In this case, there is no additional 
continuous symmetry definable and the reparametrization group SU{2)h is also contained in 50(8). Electroweak 
symmetry breaking proceeds through only one breaking pattern of the global symmetry: SO (8) into 50(7), where 
7 Goldstone bosons emerge, among them some will be absorbed by the Higgs mechanism while the remaining will 
be pseudo- Goldstone bosons. The gauge interactions respects only Gsm- Then, the breaking pattern of the gauge 
group Gsm is not defined at the tree level of the potential because one nontrivial 50(8) orbit contains both VEVs 
that are neutral and charge breaking with respect to Gsm ■ Radiative corrections of the gauge interactions should be 
incorporated to decide which model preserves U{l) em and to calculate the radiative masses of the pseudo- Goldstone 
bosons. The fate of the CP symmetry may depend on the radiative corrections of the gauge sector as well as of the 
Yukawa sector that violates CP. 

Therefore the consideration of 50(4)c as a global symmetry in NHDMs has two advantages: (i) it can preserve the 
approximate custodial symmetry 50(3)c and (ii) it avoids the presence of pseudo-Goldstone bosons and degenerate 
symmetry breaking patterns at tree level. If we insist on considering only the global 50(4)c symmetry in NHDMs, 
the immediate question towards generalization is if we can use different representations n a other than the doublet 
representation for SU{2)r. The canonical implementation uses n a — 2 for all 4> a . We immediately conclude that the 
association Y/2 — T^r is not possible for any other representation n a > 2, since we have only two hypercharges, one 
for cj) a and the opposite for </>*, common to all the doublets. For example, to have n a = 3 we need a state with zero 
hypercharge and to have n a = 4, we need two different hypercharges in modulus. Related to that, one can obviously 
embed the electroweak gauge group into a larger group that is broken at high energies to a model containing two or 
more Higgs doublets with a global accidental symmetry inherited from the high-scale symmetry 



Within the setting of Sec. Ill Dl different but equivalent implementations of SO(4)c were considered for the 2HDM 
in the literature. We reviewed them in Sec. lIIII and showed a generalization in Sec. MI CI We explicitly showed the 
basis change necessary to obtain the canonical implementation. 
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Another interesting general feature of the NHDMs implementing the global SO (4) c symmetry is the role played by 
the residual SO(N)h reparametrization freedom. For N = 2, it is just the freedom to rotate the plane perpendicular to 
kcp- For N — 3, it is the SO(3)h freedom to reparametrize t' q and its orthogonal complement. Such reparametrization 
freedom can still be used to study general SO(4)c invariant potentials. Another possible role for SO(N)h is the 
possibility to promote the whole group or some subgroup G' H to symmetries of the potential independently of the 
SO(A)c symmetry. The global symmetry in this case would be SO{A)c®G' H . Obviously, the symmetry G' H could be 
broken by EWSB. 

It is important to emphasize that the whole horizontal group SU{N)h in NHDM no longer commutes with the 
SO(A)c symmetry and some invariants that are SU(N) ^-invariant may not be S l O(4)c-invariant. That is the case 
of the invariant (r 2 — r 2 ) in 2HDM [l^, [2l[ that characterizes neutral vacuum if its null. Since it is not an SO (4) c 
invariant, its consideration alone can not distinguish among some symmetry-breaking patterns of the gauge group for 
the pattern SO(4) c -> SO(2) C - A 50(4)c-invariant measure for SO(4) c -> SO(3) c was shown in Eq. (j5Tj) . 

Another context where basis transformations are important is the connection with the Yukawa sector. Although 
all the implementations of SO(4)c symmetry are equivalent within the Higgs potential, some properties of the whole 
theory (Lagrangian), including symmetries, might depend on how the symmetry is implemented in one sector relative 
to the other. Since the global SO(4)c is explicitly broken by the gauge interactions and Yukawa interactions, the 
breaking effects depend on the implementation of the symmetry in the potential relative to the other breaking 
sectors [2^, 

Understanding how the breaking of some global symmetries occurs is also important in other contexts such as 
flavor physics. For example, Ref.[36| compares the possibility of FCNC in multi-Higgs models under the perspective 
of the two most popular suppression mechanisms, namely, Natural Flavor Conservation (NFC) and Minimal Flavor 
Violation (MFV) [371] . It is also possible to constrain 2HDM models from the requirement of mass hierarchy and 
supression of FCNC jag. 

In summary, we have shown that all implementations of the SO{A)c global symmetry in NHDMs are equivalent 
within some general conditions such as the presence of N equal representations of SU{2)r. Such results have enabled 
us to study a general S'0(4)c-invariant potential in its canonical implementation where the CP symmetry is manifest. 
We have shown that for N = 2,3 the difference between S*0 (4) ^-invariant and CP-invariant potentials is minimal 
and the knowledge of the CP symmetry immediately gives information about S l O(4)c-invariance or violation in the 
potential. We have explicitly reviewed some implementations for the 2HDM and showed a general implementation. 
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Appendix A: Custodial symmetry: SM and notation 

The Higgs potential in SM possesses a global accidental SO(A)c ~ SU(2)l®SU(2)r/Z2 symmetry due to the 
presence of one isolated doublet representation. Such symmetry group contains the electroweak gauge group Gsm = 
SU(2)l®U(1)y as a subgroup, which is the true symmetry of the Higgs Lagrangian, including the gauge interactions. 

The electroweak symmetry Gsm is then broken down to U(l) em by the Higgs VEV (<^>°). At the same time, the 
EWSB mechanism breaks SO{4)c down to the subgroup SO(3)c ~ SU(2)l+r containing U(l) em , which is the only 
nontrivial symmetry-breaking pattern possible. In the absence of U(l)y interactions [gy — > limit), the three would- 
be Goldstone bosons absorbed by the Higgs mechanism form a triplet under SO(3)c- Such symmetry protects the p 
parameter from large radiative corrections quadratic in the Higgs mass and it is called custodial symmetry 0, Q • 

Instead of distinguishing the various symmetries above by names, we denote them by their group structure and use 
the isomorphic groups interchangeably, often considering their local structure only. The groups in question are 

SU(2) L ®U(1) Y C SU{2) l ®SU{2) r /Z 2 - SO(4) c 

I EWSB I EWSB (Al) 

U(l) em c SU(2) L+R ~ 50(3) c • 

The subgroup/group relation in the first line of Eq. (|A1|) also considers U(l)y C SU(2)r. 
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Appendix B: Generators of SO(A)c 
The transformation matrix B in (|12p is given by 



B = 




(Bl) 



This relation can be read from Eq. ©. 

If we define the generators of SU{2)l+r and SU(2)l-r acting on Xa as 



the generators acting on ?y (a) are 



Or explicitly, 
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We can immediately notice that {Lj} are the generators of SO(3)c and {Kj} are the generators of the cosets 
SO(4) c /SO(3)c. 
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